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Abstract: We present a study within the fcv’-factorisation scheme on single bottom quark 
production at the LHC. In particular, we calculate the rapidity and transverse momentum 
differential distributions for single bottom quark/anti-quark production. In our setup, the 
unintegrated gluon density is obtained from the NLr BFKL Green function whereas we 
included mass effects to the Lx heavy quark jet vertex. We compare our results to the 
corresponding distributions predicted by the usual collinear factorisation scheme. The latter 
were produced with Pythia 8.1. 
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1 Introduction 

Major theoretical developments in the last three decades in small-x physics made possible 
phenomenological analyses of high energy scattering processes within the fer-factorization 
scheme [1-3] at ep (HERA) and hadron colliders (Tevatron, LHC). The Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) framework for the resummation of high center-of-mass energy 
logarithms at leading (LL) [4] and next-to-leading (NLL) [5] logarithmic accuracy is in the 
core of the majority of these analyses. 

It is very natural to wonder whether the knowledge acquired from the study of Deep 
Inelastic Scattering (DIS) processes at HERA within the BFKL formalism, mainly from 
the description of F 2 and Fi data, could be of direct use for the description of processes 
at the LHC. In principle, factorization and universality dictate the existence of a transition 
approach from ep to hadron-hadron collisions [6-8], despite the different kinematic phase 
space limits. A simple way for that to be realized and act as a proof of concept is to use 
an unintegrated gluon density from HERA fits into a phenomenological study of an LHC 
process. Recently, there were successful attempts for the detailed description of the and 
X dependence of the structure functions F2 and Fl by making use of a collinearly-improved 
BFKL equation at next-to-leading logarithmic (NLx) accuracy [9, 10].^ 

Within high energy factorization, the description of any hard process requires three in¬ 
gredients: the universal BFKL gluon Green’s function which resums high energy logarithms 
and two process dependent impact factors which describe the coupling of scattering particles 
to the gluon Green’s function. In the present case, only one impact factor (the ‘heavy quark 
impact factor’) is characterized by a hard scale i.e. the heavy quark mass and large trans¬ 
verse momentum which enables us to calculated it using perturbative QCD and collinear 
factorization. The second impact factor (the ‘proton impact factor’), which describes the 
coupling of the gluon Green’s function to the proton, is intrinsically non-perturbative and 
needs to be modeled. Combination of the gluon Green’s function and the proton impact 
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^See also the works in Refs. [11, 12]. 





factor yields then the above mentioned unintegrated gluon density. The impact factors for 
gluons and massless quarks have been calculated in Ref. [13-15], at NLx. The NLx impact 
factor for a massive quark in the initial state has been calculated in Ref. [16, 17]. 

In the last years, studies of BFKL evolution were mainly focused on processes with two 
hard scales of similar sizes in the final state to suppress any collinear-like evolution, with 
Mueller-Navelet jets [8] the best known example. Most of the studies were carried out at 
NLx accuracy [6, 18-25]. 

On the other hand there has been also considerable interest in the study of pro¬ 
cesses with one hard scale, which involve unintegrated or Transverse Momentum Dependent 
(TMD) parton density functions (PDFs). Examples of such processes at the LHC include 
forward jet [26-29] and forward Z production [30, 31]. During recent years the study of 
TMD PDFs has become a very active area of research, which find applications in various 
multi-scale processes in hadronic collisions, see Ref. [32] for a recent review. Extraction 
of TMD PDFs has in some cases been developed to very sophisticated levels, including a 
detailed discussion of experimental uncertainties, see e.g. [33]. 

In this paper, we study single bottom (or anti-bottom) quark production at the LHC. 
Bottom quark production (more accurately, bottom pair production) has received lots of 
attention in the literature [34-45] both in the collinear and the /cj’-factorization approach. 
Bottom quarks can generally be produced via gluon splitting, g ^ hb 'm. proton-proton 
collisions. Since our main purpose here is to test the unintegrated gluon density from the 
HERA fit [9, 10] and to compare it to theoretical predictions from collinear factorization 
at small x, we concentrate in the following on bottom quark production in the forward 
region of one of the protons. In this way the heavy quark - as an incoming parton - will 
be fixed at relatively large x, while the second parton - a gluon - is forced into the small-x 
region. Measurement of such a process will be possible within the LHCb experiment [57] 
and currently discussed forward updates of the ATLAS and CMS experiments. 

While the unintegrated gluon density extracted from [9, 10] does not provide a detailed 
discussion of experimental uncertainties (unlike e.g. [33]), it is the only currently available 
unintegrated gluon density which is subject to BFKL evolution at NLL accuracy, including 
a resummation of large logarithms at the level of the next-to-leading order BFKL kernel. 
In this sense the current studies present an advance over previous attempts based on LL 
accuracy. 

The article is structured as follows: in Section 2 we introduce the high energy factoriza¬ 
tion framework we will use and in Section 3 we derive the master formula for the differential 
single bottom quark cross-section. In Section 4 we present the numerical results and we 
conclude in Section 5. 

2 Forward single bottom quark production in high energy factorization 

In the following we will study for typical LHC center-mass-energies ^/s = 8 and 13 TeV the 
process 


proton(pi) -|- proton(p 2 ) bottom quark jet(A:) -|- X 


(2.1) 


where the jet rapidity is assumed to be close to the forward region of the scattering proton 
with momentum pi. We assume in the following light-like proton momenta pi and p 2 , with 
2pi ■ P 2 = s. For the above process, the bottom quark jet provides a hard scale, both 
through the bottom mass mt, and its transverse momentum kT, which allows for an an 
analysis of this process within QCD perturbation theory. Furthermore, since the scattering 
proton with momentum p 2 is separated from both the heavy quark jet and the proton with 
momentum pi by a large interval in rapidity, a description of the process within high energy 
factorization is possible. For sufficiently high kx, this process is then described at leading 
order through the partonic process Q + g ^ Q' convoluted with corresponding gluon and 
heavy quark distribution functions. Within high energy factorization, the initial heavy 
quark is always taken at large xq ~ 1, while the gluon is pushed into the small <C 1 
region, with the opposite conhguration {xq <C 1 and Xg ~ 1) suppressed by powers of 
the center-of-mass energy. For the further analysis within high energy factorization (which 
includes a resummation a large terms (aglnl/xg)” ~ 1 to all orders in 0 *), it is then 
sufficient to analyze the process 

Q{xq ■ pi) + p{p 2 ) —t bottom quark jet(A:) -|- X', (2-2) 

i.e. we study scattering of a heavy quark on a proton together with production of a heavy 
quark jet in high energy limit, see Fig. 1. The cross-section for for this process aq can be 
written as a convolution of three objects: the partonic heavy quark impact factor, the gluon 
Green’s function, which is a process independent universal quantity and the proton impact 
factor. Formally, this means that we can write for the forward bottom quark cross-section: 



Figure 1. Single bottom quark production in high energy factorization. The cross-section is 
given by a convolution of the gluon Green’s function, the proton impact factor (at the bottom of 
the diagram) and the bottom quark impact factor. A generic order heavy quark impact factor is 
depicted to the left whereas to the right the Lx impact factor is shown. 
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where Q is the hard scale related to the hnal state heavy quark momentum. We have 
introduced qi = i = 1,2 (the transverse momenta of t-channel gluons, see Fig. 1). 

In the equation (2.3) ^Q{qi,Q‘^) is the heavy quark impact factor, <hp((? 2 ,Qo) proton 
impact factor and gi, ^ 2 ) the gluon Green’s function adapted for DIS-like kinemat¬ 
ics. ^Q{qi,Q‘^) and gi, ^ 2 ) are quantities, that are calculable in perturbative QCD 

whereas ^^^(^ 2 , Qq) is an object of intrinsic non-perturbative nature and has to be modeled. 
We will use in this study the fit of Refs. [9, 10] which achieves a successful description of 
F 2 and Fl HERA data with a very simple ansatz for the proton impact factor with three 
independent parameters. 

When the two scales and Qq are similar in size, the gluon Green’s function F - 
which is obtained as the solution to the BFKL equation - can be written at leading order 
as 


( 5 , 91 , 92 ) 


2-Kqiq2j 2TTi J 2m\ql) \q1q2J w - (7) ’ 


(2.4) 


with 6is = asNc/ir and Xoil) = 2'0(1) — '(/’(t) ~ ~ t) fhe eigenvalue of the Lx BFKL 

kernel with "0(7) is the logarithmic derivative of the Euler Gamma function. The gluon 
Green’s function is universal and resums d” log” s terms to all-orders in the strong coupling. 

In our setup however, S> Qq and this expression should be written in a form 
consistent with the resummation of dg log ( 1 /x) contributions: 


-i^(s,9i,92) 


1 fdxfqfVfsY 1 

2TTql J 2m J 2T:i \ql ) \ql) ^ - asXo 


(2.5) 


In the limits 7 —)• 0,1, the zeros of the denominator of the integrand generate all-orders 
terms not compatible with DGLAP evolution [46, 47]. By taking into account the NLx 
correction to the BFKL kernel, the first of these pieces (0(a^)) is removed. Higher orders 
though however remain and are numerically important. A scheme to eliminate these spu¬ 
rious contributions was introduced in ]46] by using a modified BFKL kernel in Eq. (2.4) 
incorporating the change Xo( 7 ) ^ 2 ^( 1 ) — ' 0(7 + j) ~ 0(1 ~7 + j)- 

The NLx kernel after collinear improvements can very well be approximated by breaking 
the transcendentality of the NLx kernel and solving it pole by pole and summing up the 
different solutions. This procedure was introduced in Ref. [47] and we refer the reader there 
for further details. The NLx kernel with collinear improvements we will be using hereafter 
reads ^ 

X (7) = dsXo (7) + dsXi (7) - (7) Xo (7) + XRG{as,X, a, b). (2.6) 

with 


XRGias,X,a,b) = ds(l-h ad^) (^(7) - 0(7 - ^d^)) 
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For the NLx BFKL kernel we have: 


^ ^ c / ^ /3o 2^ ^ , ^"(7) + ^"(1 - 7) - </>(7) - </>(! - 7) 

Xi(7) = 5 x 0 ( 7 ) - —Xo(7) +-- 
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whereas the coefficients a and b read 
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To achieve a model with sensible parameters for the proton impact factor dominated by the 
infrared region, the Brodsky-Lepage-Mackenzie (BLM) optimal scale setting scheme [48] 
has been used in [9, 10] to fix the renormalization scale. ^ The BLM procedure is a way 
of absorbing the non conformal terms of the perturbative series in a redefinition of the 
coupling constant, to improve the convergence of the perturbative series. Practically, one 
needs to extract the /3o-dependent part of an observable and choose the renormalization 
scale such that this part vanishes. In the current case this leads to 


as (Wo, 7 ) = 


4Nr 


Po log (^) + ixo(7) - I + 2 (1 + |y) 


( 2 . 11 ) 


where we are using the momentum space (MOM) physical renormalization scheme based 
on a symmetric triple gluon vertex [52] with Y ~ 2.343907 and gauge parameter ^ = 3 
The modifications we need in the BFKL kernel in order to introduce this new scheme are 
as —t as iQQo,^) and Xi(7) Xi(7) ffi Eqs. (2.6), (2.8) together with the corresponding 

adjustments for the coefficients a, 6 —)• a,b which enter Eq. (2.7). The modified quantities 
read 
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^The first application of the BLM scheme was in Ref. [49-51] in the context of virtual photon-photon 
scattering. 
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where S = ^2 ^ addition, in order to access the region of small photon virtualities, 

in [9, 10], a parametrization of the running coupling introduced by Webber in Ref. [53] has 
been used. 
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(2.15) 


with A = 0.21 GeV. At low scales this modified running coupling is consistent with global 
data of power corrections to perturbative observables, while for larger values it coincides 
with the conventional perturbative running coupling constant. 

Let us add here, that in a future analysis we plan to investigate effects related to the 
choice of the renormalization scale and the choice of the parametrization of the running of 
the strong coupling (see Ref. [17] and also Refs. [54-56]). 


3 The differential cross-section with bottom mass effects included 


As already mentioned in the previous section, the non-perturbative proton impact factor 
has to be modeled. We use here the same functional form as in Refs. [9, 10]: 

which introduces three independent free parameters and has a maximum at = 6 Qq. Its 
representation in 7 space reads 

The values of the parameters Qo,5 and C were determined from a fit to combined HERA 
data. When the leading order photon impact factor was used the obtained values were 
Qq = 0.28 GeV, 6 = 8.4 and C = 1.50 whereas in the case of the kinematically improved 
photon impact factor the last two change to C = 2.35 and h = 6.5 with the number of 
flavors fixed to n/ =4. We use both sets of values for C and 6 in our numerical study later. 
Combining the BFKL Green’s function for DIS kinematics and the proton impact factor, 
we obtain the following expression for an unintegrated gluon density within our setup 

G{x,qi)= [ ix,qi,q 2 )<^p{q 2 ,Ql) ■ (3.3) 

J Q 2 


To obtain the complete NLx BFKL Green’s function we need to add to Eq. (2.4) apart from 
the NLx correction to the BFKL eigenvalue, non-exponentiating NLx /3o terms. Following 
the treatment of Ref. [10] one obtains 


G{x,qi,Q) = 


qt J-c 


du 

2tt‘^ 
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In the DIS analysis has been identified with the virtuality of the photon. In the present 
study we use instead the transverse momentum of the bottom quark, Q = kx, with kx = 
the modulus of the transverse momentum of the heavy quark. The other obvious 
choice for Q = \Jk‘^ + rri^ causes only small differences in the results. 

Once we have the formal definition of the gluon density given by Eqs. (3.3) and (3.4), 
we can have an alternative view at QP’Scattering depicted in Fig. 1. In particular, we 
may consider it as a convolution of the bottom quark impact factor with the gluon density 
as shown in Fig. 2. To obtain the complete forward heavy quark cross-section, we further 



Figure 2. The same process as in Fig. 1 here presented as a convolution of the bottom quark 
impact factor and the gluon density. 

require the bottom quark jet vertex which at leading order coincides with the massless quark 
impact factor modulo Dirac Delta functions to ensure momentum conservation. The heavy 
quark jet vertex at Lx depends therefore only implicitly on the quark mass through the final 
state phase-space integration. We should stress here, that a proper NLx study with mass 
effects properly introduced would be the desirable goal. Currently only the inclusive heavy 
quark impact factor, which describes the process Q + g* ^ X", is available at NLx accuracy 
[16, 17]. While at leading order inclusive and jet impact factor coincide, a NLx description 
of the process Q + g* ^ bottom quark jet -|- X"' will depend explicitly on the details of the 
employed jet algorithm and hence differs from the corresponding inclusive result. In the 
current study we therefore restrict ourselves to the Lx jet vertex. The leading order impact 
factor is then obtained from the squared amplitude of the subprocess Q+g* —)• Q', integrated 
over the one-particle invariant phase space The momenta of the incoming quark and 

gluon can be expressed in Sudakov variables as pQ = xgpi and pg = XgP 2 + Qi respectively. 
Energy-momentum conservation identifies then the gluon transverse momentum with the 
transverse momentum of the final state bottom quark qi = k\ the one-particle phase space 
of the heavy quarks reads = 2'K5{xqXgS — rn^). After a bit of Algebra, we obtain 
for the Q + g* ^ Q' partonic cross-section: 

- r/ 1 2 2\ (y.g^'K 

a = aod{xQXgS - k - m^) , ao = ^ ■ 
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The total pp ^ Q + X cross-section then reads 


^pp^Q+X 



[fQ{xQ,IJ^f) + fQ{xQ,fJ^f)] G{xg, kr, Q) . 


(3.5) 


with fi, i = Q,Q the collinear (anti-) bottom quark distribution and pf the collinear 
factorization scale. Fixing Xg = {k?‘+m‘l)/xQS and introducing the rapidity of the produced 
bottom quark rj = ^In^, the total cross-section is recast into 


’Xpp^Q+X 



cPk 

TT 

cPk 

TT 


- XQ [fQixQ,pf) + fQ{xQ,pf)]) G{xg, kr, Q) 

XQXgS 

,2 1 ° 2 [fQixQ, fj-f) + fQ{xQ,fJ-f)] G{xg, kr, Q) ■ 

Krp Tlfl^ 


(3.6) 


After integrating over the azimuthal angle of fc, the pp ^ Q' X double differential cross- 
section hnally reads 


dcTpp^Q+X 
dr] dkx 


2kj' ■ (To 

k^ + ml 


XQ [fQ{xQ,pf) + fQ{xQ,pf)] G{xg,k‘^,Q) . 


(3.7) 


with 




(3.8) 


Leaving aside for the time being the dependence on the collinear bottom quark distribution 
function xq fi{xQ), i = Q,Q we will have for the Qp —)• Q' cross-section in zz-space 


daqp^Q' _ (JQ-C 7 ^ T {6 -l-iu) / ^ 

drjdkT kT-{k‘^ + ml) J vr^ r((5) VQo/ 


X 


|l log 







(3.9) 


Fixing the factorization scale of the bottom quark PDF to pf = y k'^ + ixil, the double 
differential cross-section for single bottom quark production in proton-proton collisions 
reads 

, (3.10) 

which we use to produce all of our numerical results in the next section. 


dcr, 


PP^Q'+X 

dr] dkx 


dcTQp^Q' 

dpdkx 


XQ 


fQ[xQ,\jk‘^ + ml \ + /q( + m 


4 Numerical results 

In this section we present our predictions for the differential cross-section, p- and k^- 
distribution, for single bottom quark/anti-quark production at the LHC. The analysis does 
not make distinction between the bottom quark and anti-quark. 























In Fig. 3, we show the range of maximal/minimal Xg for each value of kx- The light 
blue area corresponds to the Bjorken x and Q ranges of the data, that were used for the 
F 2 and Fl fit in [9, 10]. The light orange area corresponds to the Xg and kx ranges in our 
calculation of the bottom quark cross-section when the rapidity takes values between 1 and 
5 and kx is constrained to 4 GeV < kx < 100 GeV. Lastly, the light red area corresponds to 
the X and kx ranges when the rapidity takes values between 3 and 5 and kx is constrained 
to 4 GeV < kx < 50 GeV. This kinematical range is not covered by the general purpose 
detectors ATLAS and GMS, but is accessible by the LHCb detector [57] designed for these 
kind of measurements. 

It is more than evident, that the kinematic region we are covering for bottom quark 
production does not overlap at all with the kinematic region in which the unintegrated gluon 
density was obtained. In particular, the values of the unintegrated gluon density tested in 
our setup are not directly constrained by HERA data. They are rather calculated through 
evolving the results of the HERA fit towards both smaller x values and larger kx using 
our collinear improved solution to the BFKL equation. Nevertheless, as we shall see, the 
unintegrated gluon density based on that specific model, gives results very close to Pythia 
8.1. It has to be stressed though, that the whole approach carries uncertainties which at 
present cannot be quantified and one has to be cautious not to interpret these high energy 
factorization results as the final word within the BFKL approach. 

In all the next figures, predictions with 6 = 8.4 and C = 1.5 are plotted with solid 
red lines, the results with 6 = 6.5 and C = 2.35 are plotted with dashed red lines and the 
results by Pythia 8.1 are plotted with purple solid lines. We have used the MSTW 2008 
NLO parton density functions [58] throughout the entire section. The bottom quark mass 
was set to mb = 4.7 GeV for the high energy factorization result whereas in PYTHIA 8.1 
the program default value was used. 




Figure 3. Light blue area: kinematic region corresponding to the original fit of the proton impact 
factor. Light orange area: phase space for the kinematic cuts (1 < 77 < 5, 4 GeV< kx < 100 GeV). 
Light red area: phase space for the kinematic cuts (3 < 77 < 5, 4 GeV< kx < 50 GeV). 

In Figs. 4 and 5 we present results for y/s = 8 TeV and in Figs. 6 and 7 for y/s = 13 
TeV. In Figs. 4 and 6 we integrate the differential cross-section (3.10) over the rapidity of 
the quark in the range 1 < 77 < 5 and over kx in the range 4 GeV < kx < 100 GeV while in 
Figs. 5 and 7 we integrate over the rapidity of the quark in the range 3 < r/ < 5 and over the 
transverse momentum in the range 4 GeV < kx < 50 GeV. The differential distributions are 
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Figure 4. Collision energy y/s = 8 TeV. Left: 77 -distribution after integrating over kr in the range 
4 GeV < kr < 100 GeV. Right: fc^-distribution after integrating over 77 in the range 1 < 77 < 5. 
Both distributions are normalized by the integrated cross-section over 77 and kx in the ranges 
4 GeV < fcr < 100 GeV and 1 < 77 < 5. 




Figure 5. Collision energy ^ = 8 TeV. Left: 77 -distribution after integrating over kr in the range 
4 GeV < kx < 50 GeV. Right: /c-r-distribution after integrating over 77 in the range 3 < 77 < 5. 
Both distributions are normalized by the integrated cross-section over 77 and kx in the ranges 
4 GeV < kx < 50 GeV and 3 < 77 < 5. 




Figure 6. Collision energy y/s = 13 TeV. Left: 77 -distribution after integrating over kx in the range 
4 GeV < kx < 100 GeV. Right: /cr-distribution after integrating over 77 in the range 1 < 77 < 5. 
Both distributions are normalized by the integrated cross-section over 77 and kx in the ranges 
4 GeV <kx< 100 GeV and 1 < 77 < 5. 


normalized by the integrated cross-section over the corresponding rapidity and kx ranges. 

A first observation is that when we compare the kx distributions calculated in high 
energy factorization and by Pythia 8.1 we see the former to be smaller than the latter due 
to the difference in the shape at small kx- The shape of the kx distribution though is very 
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Figure 7. Collision energy ^/s = 13 TeV. Left: r^-distribution after integrating over kr in the range 
4 GeV < /ct < 50 GeV. Right: /cT-distribution after integrating over 77 in the range 3 < 77 < 5. 
Both distributions are normalized by the integrated cross-section over rj and kx in the ranges 

4 GeV < fcr < 50 GeV and 3 < 77 < 5. 

similar at large kx in both approaches. 

The main finding comes forward when we focus on the rapidity distributions. The high 
energy factorization result seems to be somehow larger than the Pythia 8.1 estimate for 
small rapidities whereas for larger 77 it drops faster than the Pythia 8.1 result and at the 
high end of the rapidity it lies below it. The same trend is followed for both center-of-mass 
energies and for both kx integration ranges (4 GeV < /cr < 50 GeV and 4 GeV < kx < 
100 GeV). 

Finally, let us note, that a change of values for the parameters 5 and C from 5 = 8.4 
and C = 1.5 to 5 = 6.5 and C = 2.35 results to a slightly smaller NLx cross-section. 

5 Conclusions and Outlook 

We have presented a study of the rapidity and kx (transversal momentum) differential 
distributions for single bottom quark production at the LHG calculated both in high energy 
factorization and by the Monte Garlo program PYTHIA 8.1. and for ^/s = 8 and 13 TeV. 
Within the former framework, we have used a model for the proton impact factor, the NLx 
BFKL gluon Green’s function and the Lx heavy quark jet vertex with bottom mass effects 
included. 

The main result of our study concerns the rapidity distributions. The high energy 
factorization estimate is for small rapidities larger than the PYTHIA 8.1 result but as the 
rapidity approaches some middle range value, the fall becomes steeper and the estimate 
gets smaller than the PYTHIA 8.1 estimate. The /c-r-distributions have very similar shapes 
but the Pythia 8.1 result is almost always larger. 

The calculation presented in this article suggests, within its limitations, namely, that 
it is only a partial NLx calculation and that the unintegrated gluon density is probed in 
a region not covered by the F 2 and Fx fit, that single bottom quark production might 
be used as an experimental probe of the /c^-factorization scheme and the validity of high 
energy factorization for LHC processes. It also shows, that our initial assumption, that the 
unintegrated gluon density from HERA would not fail at the LHC, was justified. We plan to 
extend our study toward obtaining more exclusive information with regard to the final state 
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by using the BFKL Monte Carlo code BFKLex [59, 60]. The code is an implementation of 
an iterative solution to the NLx BFKL equation and has already been used in a number of 
projects [61-63]. We expect, that once the NLa; massive quark jet vertex is ready, we will 
able to present a more refined study of the process and a detailed analysis on the BFKL 
predictions for the single bottom quark cross section. 
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